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ABSTRACT 

The cores of luminous B and A-type (BA) supcrgiant stars are the seeds of later core collapse super- 
novae. Thus, constraining the near-core conditions in this class of stars can place tighter constraints 
on the size, mass and chemical composition of supernova remnants. Asteroseismology of these massive 
stars is one possible approach into such investigations. Recently, Moravveji ct al. (2012, hereafter Pa- 
per I) extracted 19 significant frequencies from a 6-year radial velocity monitoring or Rigel (/3 Ori, B8 
la). The periods they determined broadly range from 1.22 to 74.74 days. Based on our differentially 
rotating stellar structure and evolution model, Rigel, at its current evolutionary state, is undergoing 
core He burning and shell H burning. Linear fully non-adiabatic non-radial stability analyses result in 
the excitation of a dense spectrum of non-radial gravity-dominated mixed modes. The fundamental 
radial mode {£ = 0) and its overtones are all stable. When the hydrogen burning shell is located 
even partially in the radiative zone, a favorable condition for dcstabilization of g-modcs through the 
so-called e-mechanism becomes viable. Only those g-modcs that have high relative amplitudes in the 
hydrogen burning (radiative) zone can survive the strong radiative damping. From the entire observed 
range of variability periods of Rigel (found in Paper I), and based on our model, only those modes 
with periods ranging between 21 to 127 days can be theoretically explained by the e-mcchanism. The 
origin of the short-period variations (found in Paper I) still remain unexplained. Because Rigel is 
similar to other massive BA supergiants, we believe that the e-mechanism may be able to explain the 
long-period variations in a Gygni class of pulsating stars. 



1. INTRODUCTION 

After the core hydrogen depletion, massive stars (M > 
15Mq) enter the BA supcrgiant phase of stellar evolu- 
tion. During this epoch of evolution, they exhibit peri- 
odic and/or quasi-periodic low amplitude flux microvari- 
ability in addition to variations i n radial velocity and 
equivalent wi dth of spectral lines fWaclkcns et al.lll998l : 
iKaufer et all [1996, 19 97; Bres olin et al . 200j). Spec- 
troscopic analysis of iLefever et al] (j2007D clarified that 
blue supergiant s (BSGs) pulsate in the gravity mode . 
Recent studi es (IDegroote et all 120101 : INoels et al.l 120101 : 
IMiglio et al.ll200l " show that asteroseismology of mas- 
sive BSG stars can probe the extent of convective core 
overshooting, based on the observed and predicted regu- 
larity in period spacings. 

On the other hand, there are various theoretical ex- 
planations for the obser ved oscillat i ons in BSGs. The 
photometcric study of iSaio et all ()2006l ) on the B2 
Ib/II BSG HD 163899 (which is most likely a less- 
evolved less-massive analogue of Rigel) showed that the 
K-mechanism can excite a rich spectrum of pressure (p) 
and gravity (g) modes in the post-terminal-age-main- 
sequence stars. More importantly, they showed that only 
those modes that arrive with an appropriate phase to 
the base of intermediate convective zone (ICZ) can be 
reflected back to the surface, and hence be observed. 
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iGautschvl (|2009| ) searched for the origin of long-period 
variabilities in the prototype a Cyg. Interestingly his 
Figure 5 shows a gap for stellar models with 3.95 < 
log Toft ^ 4.15 where no instability is predicted; Rigel 
lies in this gap. According to Sajo (201 1.1. Rigel should 
be unstable against non-radial convective g~ modes. Ra- 
dial and non-radial strange mode is proposed as another 
mechanism to induce instability and interplay with mass 
loss efficiencies in those massive stars with log(_L/ M) > 
4 log( /Mp, ) ( Gla t_zcl & Kiriakidis 1993; Sai o et al. 
1998t iGlatzel et ai.l ll999l: iDziembowski &: Slawinska 



2005l: IAerts et al.ll2010HSaidl2011l ). This requirement is 
also marginally fulfilled by Rigel. iGodart et all ()2009l ) 

investigated the destructive impact of core overshoot- 
ing and mass loss during MS on the extent of ICZ, and 
showed that models with wider ICZs are more likely to 
destabilize stellar oscillations. 

Once the MS evolution of massive stars ends due to the 
H depletion in the core, the CNO fusion which takes place 
in the hydrogen burning shell (hereafter HBSh) is the 
main source of energy. Further contraction and heating 
of the helium-rich core initiates core He burning (CHcB) 
when the star is a late B supcrgiant. While the triple- 
alpha (3q;) and CNO reactions have mild p— dependence, 
they are both sensitive functions of temperature. In the 
zeroth order, the energy generation rate can be expressed 
as Cnuc pT"^ where 



(1) 
(2) 
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TABLE 1 

Observed Rigel properties and MESA Input Models. 
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Regarding the 3a reaction, Ep = 2, and €t = 40 to 19 for 
Tg = 1 to 2 in addition to an extra dependence. Sim- 
ilarly, for the CNO burning network ep = 1, and £7^ = 13 
to 17 at Tg = 50 to 25. Tg and Tg denote temperature in 
million and hundred million Kelvin degrees, respectively, 
and Y is the fractional He abundance. For more de - 
tailed discussion see Clayton (1968), Wei ss et al.l (|200l . 
IMaederi (|2009t) and fWiescher et al.l pOlOj. 

Several authors have shown that durin g the 
evolu ti on of low mass stars ([Noels et al.l 11974 
Il976t iSonoi fc Shibahashl [Ml and massive stars 
([Shibahashl &: Osakil 1 19761) . there occurs conditions 
where the Cnue can inject enough energy into the 
mechanical energy required for the propagation of 
vibrational modes, and overcome strong radiative 
damping close to the stellar cor e. This is the so called 
e-mechanism ([Unno et al.lll989l ). Th is mechanism can 
be operative duri ng shell H burning ([ Kawalerl 1 19881 ) or 
shell He burning (iKawaler et al.|[T986() phase to success- 
fully explain the observed gravity mode instabilities in 
hot white dwarf stars. The radial fundamental mode 
of very low mass stars and brown dwarfs can be also 
unstable during the central deuterium burning phas e 
([Palla fc Baraffi [2001 iRodriguez-Lopez et al.l [20Tll). 
Most recentlv IMiller Bertolami et al.l ([201 ll ) applied the 
same mechanism to a subdwarf B star prior to CHcB 
phase. In accord with these, the favorable conditions 
for excitation of internal g-modes can also take place 
in post-MS massive stars, and lead to pulsational 
instabilities. 

We find Rigel to be an excellent testbed to investi- 
gate the existence of the e-mechanism to excite oscilla- 
tory modes and to study the possible asteroseismic po- 
tential of slowl y pulsating blue superg iant stars. In a 
previous work, iMoravveii et al.l ([2012L hereafter Paper 
I) presented the MOST photometry and radial veloc- 
ity variations of Rigel, and tabulated 19 significant pul- 
sation frequencies from radial velocity variations. The 
periods of observed modes range from 1.21 to 74.74 
days. Table [T] summarizes the observed physical prop- 
erties of Rigel from the literature which are an invalu- 
able input for the equilibrium modeling of the star. Yg, 
(N/C)s, and (N/0)s arc surface abundance ratios of He 
and CNO products, respectively. The column marked 
with V shows the measured v sin i ([Przvbilla et al.l[20Tol : 
iSimon-Diaz et afl2010( ) . We explain the properties of the 
model we found to describe Rigel in Section [21 The pul- 
sation analysis is discussed in Section [3l and in Section [4l 
we summarize the conditions for the occurrence of the e- 
mechanism. The purpose of this paper is to highlight the 
importance of the e-mechanism as a possible destabiliz- 



ing driving force agains pulsations in BSC core-collapse 
progenitors ([Murphv et a l. 2004). 

2. MODELING RIGEL 
2.1. Input MESA Model 

A model that can represent the current evolution- 
ary stage of Rigel on the Hertzsprung-Russell diagram 
(HRD) has been generated with the aid of MESA stel- 
lar structure and evolution code (version 3723). A full 
descri ption of this instrument is given by iPaxton et alJ 
([20 lit ). A non-magnetic differentially rotating model 
with the initial mass Mzams = 19-M^o and i nitial so- 
lar composition as in iGrevesse fc Noelsl ([19931 ) at solar 
mctallicity Z = 0.02 were generated with the initial 
rotation velocity Vi = 200 km to help mixing He 
and CNO products to the surface and reproduce the 
observed values. Surface abundance yield s up to the 
re-co re-collapse stage were compared with iHeger et alJ 
2005[ ) to ensure for their agreement. Note that our 
choice of initia l abundances complies with those given 
in Tables 1 of IPrzvbilla et al.l ([2010( ). The surface of 
the model is located at the optical depth r = 0.66 7 
using Eddington T — r integrati on scheme lE ddingta 
Il926f ). Mixing Length theory of iBohm-Vitensq ()195 
is used to treat the convective mixing. The borders of 
convcctive zones are found using the Schwarzschild cri- 
teria (Viad = Vad)- Thc modcl was evolved with the 
above mentioned initial values until thc post-MS phase 
and stopped when Toft = 12 058 K, close to the recent 
measurement. The final radius of the model proved to 
have negligible dependence on the choice of the mixing 
length parameter (which is defined as the ratio of the 
length over which convective cells loose their identity to 
the local pressure scale height c^mlt = L/Hp). A value 
of aMLT = 1-6 was adopted for the model. The over- 
shooting from the top and bottom of convective zones 
were sup pressed. The me talicit y-dependent mass loss 
recipe of lVink et all ([200l and iGlebbeek et al.l ([20091 ) 
with 77 = 0.8 as in lMaeder fc Mevnetl ( 200ll ) is used dur- 
ing the evolution across thc HRD in all models. Calcula- 
ti qn of Brunt- Vaisala fre quency is based on Equation 14 
in iBrassard et al.l ([1991|) . To reach the surface yiel ds of 
He and CNO as measured bv IPrzvbilla et al.l ([2010D . we 
enhanced the efficiency of rotational mixing by a factor 
of two ( fc = 2 in Eq. 53 in IHeger et aI1l2000( ). 

Table [Tj lists the modeled physical parameters of Rigel 
to compare with observations. Not all calculated models 
meet their measured counterparts. It is neither easy to 
find a model that simultaneously matches all observed 
quantities and nor our goal to present the "best" model 
of Rigel. Yet, the adopted model closely represents the 
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Fig. 1. — Internal structure of the input model. On the left, the profile of total luminosity (short dashed line), radiative luminosity (solid 
line) and opacity (dotted line) along with distribution of mass (long dashed line) and radius (dashed-dot line) are plotted against logT. 
On the top right panel is the profile of enuc (solid line), gravitational energy (dotted line), neutrino energy loss (dashcd-dotcd line) and 
the net energy (dashed line). The right middle panel shows the abundance of He and CNO cycle elements. The bottom right panel shows 
the temperature (Eq. [l} and density (Eq. O derivatives of enuc, respectively. The vertical lines mark the extent of HBSh, and show the 
limits of tj- which is favorable in this zone for e-mechanism. 



current observed and physical properties of Rigel, and 
furthermore yields consistent seismic results. 

2.2. Internal Structure of the Adopted Model 

To comprehend better the complex internal structure 
of our model, Figure [T] shows profiles of several thermo- 
dynamic quantities with logT on the abscissa. On the 
left panel; the solid line shows the profile of radiative lu- 
minosity which starts to grow from the core due to the 
CHcB followed by a steep increase around logT = 7.7 
arising from HBSh; this profile deviates from total lu- 
minosity (short dashed line) wherever convection assists 
radiation in energy transport. Therefore, the dips in 
radiative luminosity profile mark the convective zones. 
Most relevant to our study are three radiative zones and 
three convective zones. They have profound impacts on 
the energy balance of pulsation modes. The ICZ is sit- 
uated just above the HBSh at 7.1 < logT < 7.6; the 
Iron- and He-bumps of opacity arc around logT ~ 5.2 
and logT = 4.6, respectively. There are two radiative 
zones below and above the ICZ, where the former hosts 
the peak of energy generated from H burning, and the 
latter lies above the ICZ. There is another radiative zone 
surrounded by the two opacity bumps. The long-dashed 
and dash-dotted lines signify the spatial mass distribu- 
tion of the model. Accordingly, around 99% of the mass 
is inside 25% of the star's radius. 

Right panels of Figure [1] show profiles of various quan- 
tities related to the energy generation rate. On the top 
panel, the solid line is the run of energy generation rate 
enuc in logarithmic scale. The CHcB peak in the core and 
HBSh peak at logT = 7.7 are evident. The dash-dotted 
line shows the neutrino cooling rate e^, which is mainly 



an order of magnitude smaller than enuc except below 
HBSh. There is an extra source (sink) of energy depend- 
ing on whether a zone is contracting (expanding). This is 
designated by eg ~ —TdS/dt (dotted line) where T and 
dS/dt denote temperature and entropy time derivative, 
respectively. Note that for the layers undergoing expan- 
sion, Cg is negative. The radiative zone above the ICZ is 
where expansion is dominant. Finally, the profile of the 
generated net energy Cnet is shown by the dashed line. 
Note that in the non-burning envelope, expansion and 
contraction are sources and sinks of energy. Two vertical 
dotted lines show the extent in the model where energy 
is purely generated by CNO cycle. A fascinating feature 
in this model is that this shell is situated partially in a 
radiative zone which is central to our study. 

The right middle panel of Figure [T] shows the runs 
of ^He, i^c, and ^^O abundances. The ratios of 
i2c/i4n= 0.015 and i6q/ i%= q gp comparable 
with those from Table 25.3 of lMaeded ()2009D in the HBSh 
zone. 

The right bottom panel of Figure [1] shows the logarith- 
mic derivative of energy generation rate e^ (solid line, 
Eq. [1]) and ep (dotted line, Eq. [2]) for our model, re- 
spectively. Consistently, er varies from ~30 to 40 in the 
Helium rich core, and steeply declines to ~11.5 to 14.0 
in the HBSh where CNO is contributing to a large frac- 
tion of Rigel's luminosity. This latter range of er is one 
of the favorable conditions for e-mechanism to operate 
(Section 13. 4p . A caution must be given to a peak in er 
at non-burning envelope around logT = 5.3 where Cnuc is 
vanishingly small. This arises from careful numerical dif- 
ferencing followed up in MESA and will have no effect on 
our modal stability analysis since the combination enuc 
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Fig. 2. — Propagation diagram for the model described in Table 
[T] The Brunt- Vaisala A'^ (solid line) and the Lamb frequency 5^—1 
(dashed line) are in logarithmic scale. The radiative (convcctive) 
zones are where A'^^ > (Af^ < 0). The profile of enuc is in an 
arbitrary unit only to show that the peak of enuc at HBSh resides 
in the innermost radiative zone. Four selected frequencies from 
(Table 2 in) Paper I are also shown as horizontal dotted lines to 
mark the propagation cavity of each mode. 

times eT play the major roles. 

2.3. Propagation Cavities 
The Brunt- Vaisala frequency N"^ which is defined as 



(jBrassard et al.lll99H) 



1 d\x\P 9 In p 



Fi 91nr Slnr 



(3) 



shows the extent and boundaries of convective zones 
(7V^ < where g-modes arc evanescent) and the radiative 
zones (7V^ > where g-modes are oscillatory). The com- 
plexity of the internal structure of our model is indicated 
in Figure [5] which shows the run of the square root of N"^ 



(solid line) and Lamb frequency 5| 



a function of temperature. Here, Cs denotes the sound 
speed in the stellar plasma. The profile of e^uc is desig- 
nated by dot-dashed line in arbitrary units. 

The remarkable feature in the model is that the peak 
and hotter part of HBSh reside in the innermost radiative 
zone and extend slightly to the base of ICZ. The coinci- 
dence of this condition is one of the prerequisites for the 
e- mechanism to operate (Section 13. 4p . This, however, 
is not true during the entire evolution of a massive star. 
Once the central H burning is ceased, the helium rich core 
starts to contract. As a result the core temperature Tc 
increases, and the opacity k near the core drops. In addi- 
tion, the non-burning core contributes very poorly to the 
radiative luminosity L^. Therefore, the radiative tem- 
perature gradient Vrad KPLr/mT^ becomes smaller 
than Vad = {d\nT / d\n-P)s mainly because of a dras- 
tic decrease in Lr/m. Thus, the core becomes radiative. 
This marks the starting point of the post-MS evolution- 
ary phase. The HBSh, which at this evolutionary stage, 
is the main source of energy, moves away from the core 
and lies partly in the ICZ and partly in the innermost 
radiative zone. 

3. SEISMIC ANALYSIS 



The model discussed in Section [5] is the i nput to the 
Granada oscillation Code (h ereafter GraCo: iMova et al.l 
[200llMova fc Garridol[2008l) . The radial and non-radial 
oscillation spectra are calculated in both the adiabatic 
and non-adiabatic frames sol ving numerically th e pertur- 
bative equations described in lUnno et al.l ()1989[ ). The set 
of adiabatic equations provide a first list of eigenfrequen- 
cies and eigenfunctions. In the non-adiabatic approx- 
imation, the eigenfrequencies are corrected (slightly for 
non-adiabatic effects), and an additional set of eigenfunc- 
tions related with the energy interchanges are obtained. 
These eigenfunctions make it possible to calculate the- 
oretically the work done by the pulsational mode to its 
surroundings, during a complete period. In addition, we 
can obtain the relative flux variations and the phase-lag 
between the flux variations and the corresponding varia- 
tions in radius. These last two non-adiabatically derived 
quantities are not used in the present study, but could 
be useful once multi-wavelength photometry of Rigel is 
carried out. 

In the GraCo, the non-adiabatic calculations can be 
done to include the convection - pulsation interaction 
using the Time De pendent Convection (TDC) theory 
(jPupret et al.l 120051 ) . This theory improves the "frozen 
convection" (FC) approximation usually implemented in 
most of the other codes. But this is not necessary in 
this study, since near the core of stars like Rigel, where 
the e-mechanism is located, convection is highly devel- 
oped. As a result, convection can be regarded decou- 
pled from pulsation. In this case the FC is an appro- 
priate approximation. On the other hand, the adiabatic 
solutions of GraCo have been used as reference for the 
EST A works (Evolution and Asteroseismic Tools Activ- 
ities, iLebreton et"all[200l iMova et al.l[2008l ). 

To comply with our list of observed frequencies pre- 
sented in Paper I, we restrict our search for frequen- 
cies from less than 10~^ to 1.3 d~^ and for radial and 
non-radial modes up to ^ ~ 3. Below the lower fre- 
quency limit of 10~^ d~^, convergence is very difficult 
to achie ve. In this respect, o ur strategy is different from 
that of iGodart et al.l ()200 9f) where they employ quasi- 
adiabatic approach for stability analysis of the radia- 
tive core and full non-adiabatic solution for the enve- 
lope. This approach is optimized to extract a discrete set 
of excited modes reflected from the edges of the ICZ in 
those models where the K-mechanism is responsible for 
pulsat ional instabilities. We note that while .Gau tschyl 
(|2009| ) did not find any i nstab il ities f or models of Rigel 
(3.95 < logTeff < 4.15), (I20T1I) predicted that a 
Cyg type pulsating stars ought be unstable against con- 
vective gravity modes. In the following, we explain the 
asteroseismic results for different degrees I. 

3.1. Energy Balance 

The temporal variation of perturbed quantities in 
GraCo are treated as exp(— icrt) where and ct/ are the 
real and imaginary parts of the complex eigenfrequency 
c = f^R + *cr/, respectively. The imaginary part is essen- 
tially obtained as a solution of the system of differential 
equations, however, it is also related to the so-called cu- 
mulative work W . 
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Fig. 3. — Loft. Variation of growth rate r] with mode frequency for < £ < 3. All radial p-modes (empty circles) are stable. 
The low- frequency dipole (filled circles), quadrupole (filled squares) and octupole (filled diamonds) modes are found excited by the e- 
mechanism (rj > 0). The period of these excited modes vary from about 127 days down to 21 days. There is a peak in all stable non-radial 
modes around 0.26 d~^ where a tendency towards overstability appears regardless of degree. The vertical dotted lines mark the observed 
frequencies from Paper I. Right. Similar to the left, but with the periods on the abscissa. With the increase in degree £, the transition 
from overstable to stable modes shifts towards shorter periods. 

of enuc by er (see Eq. [6|). 

3.2. Frequency Domain of Excited Gravity Modes 
In evolved stars with alternating radiative and convec- 



where Ek is the kinetic energy. According to lUnno et alj 
(|1989f ) (their Section 26), the work integral accumulates 
the contribution from perturbation to radiative flux Wf, 
convective flux Wc , and energy generation rate , say 



W ^Wf + 

The net contribution of e^uc 
mode is 



Wc + W„ (5) 
to the work integral of each 



CTi? Jo 
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where et and ep are as in Eq. [T] and [21 — 1 = 
{d\nT / d\n p)s is an adiabatic exponent and 5T is the 
Lagrangian temperature perturbation. As the bottom 
right panel in Figure [1] shows, et and ep/(F3 — 1) are 
both positive quantities; as a result, the integrand in Eq. 
[6] is greater than zero and always has a tendency towards 
mode excita.tion g iven hig h enough values for ey >11 
(|Noels et al.lll97¥ fl976: S onoi fc Shibahashi|[20Tlh . 



We define the e— folding time = aj^ , and the growth 
rate 77 = W/\W\ where \W\ is an absolute value of the 
work. With these, growth rate lies in the range — 1.0 < 
77 < 1.0. Therefore, the stability of non-radial modes 
is determined by the sign of growth rate where positive 
(negative) values denote overstability (stability). 

Since enuc is practically (but not necessarily numeri- 
cally) zero in the envelope of our model, the rise in the 
profile of Et in (bottom right panel of Figure [1} has no 
impact on the evaluation of due to the multiplication 



tive zones and large density contrast between the core 
and the envelope, pressure (p) and gravity (g) modes 
appear in mixed character, where for each, there is a 
contribution from the other. Consequently, the modes 
are identified as p- or g-modes according to the more 
domi nant effect one might have over the other (jOsakil 
[T975I) . In our case, adiabatic eigenfunctions are high- 
order gravity-dominated mixed modes with large num- 
ber of nodes in both pressure and gravity cavities. The 
huge number of nodes in the e-mechanism trapping zone 
is met by the grid of mesh points whi ch is fuhy in 
agreement with the asy mptotic approach ()Tassoullll980l : 
ISmevers fc MoTOl[2007li . 

To study the effect of energy balance. Figure [3] shows 
the variation of growth rate 77 as a function of mode fre- 
quency (left panels) and period (right panels) for dif- 
ferent spherical harmonic degrees i. Because models of 
evolved stars have highly non-adiabatic envelopes, the 
numerical convergence is not straightforward in all cases, 
so wc do not assume that all eigensolutions are covered; 
however, the solution space is well sampled so that there 
is a clear trend in stability behavior of non-radial modes 
in Figure [3] 

For non-radial modes (panels b, c and d in Figure [3]) 
and at the lower frequency range, the dense spectrum of 
eigenmodes are excited by the e-mechanism (discussed in 
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Fig. 4. — An example of two g-modes, gov (left) and ggt (right). Concise mode information is also included: £, u, P and r] are degree, 
frequency in d~^, period in days and growth rate, respectively. The modulus of the radial displacements ^rC*") in panels (a) are normalized 
to unity at the surface. Though the e-mechanism in both modes is the main driver of g-mode pulsation, the latter is stabilized by radiative 
damping above ICZ. The excited mode gov has a large amplitude in the radiative HBSh and above ICZ so that radiative damping is 
overcome. The middle panels (b) show the work integral W (solid black lino), work derivative dW/dlnr (dotted red line) and square of 
entropy perturbation (SS)'^ (dashed blue line) all normalized to unity. The bottom panels (c) are zoomed-in views of the same quantities 
as in panels (b) in the HI3Sh. 



more detail in Section [3^ and the overstability declines 
steeply by frequency with an abrupt transition from un- 
stable to stable modes. At the high frequency end, all 
modes are strongly damped rj w —1.0. In between these 
two extremes there exist a broad peak in the stable non- 
radial modes occurring at --0.26 d"^ This shows that 
there is a tendency towards instabilities at higher fre- 
quencies whose driving force is different from the low- 
frequency modes. As it is discussed in more detail in 
Section 13.51 the inefBciency of the ft-mechanism oper- 
ating around the He bump of opacity is responsible for 
their stability. Very clearly, there are many observed fre- 
quencies (marked by vertical dotted lines) that have no 
predicted unstable counterparts. 

Panels (e) to (h) in Figure [3] present the profile of the 
growth rate with the corresponding period of the modes. 
The excited modes are more distinct in this represen- 
tation. It is clear that four observed modes lie in the 
predicted domain of excited g-modes destabilized by the 
e-mechanism. 

An intriguing question would be the dependence of 
the frequency domain of overstable modes on the in- 
put physics of the equilibrium model. In this way, the 
low-frequency stable modes could be helped into oversta- 
bility. Such physical ingredients could be the increased 
opacity at Carbon bump (logT w 6.3), efficiency of rota- 
tionally induced mixing, overshooting from the borders 



of the ICZ and the mass loss rate. We intend to investi- 
gate this in the upcoming studies. 

3.3. Radial Modes 

The Fundamental radial mode (pi having a period of 
8.02 d) and it's overtones in the prescribed scanning 
frequency range are very stable. This stability takes 
place around the Helium partial ionization zone close to 
the surface o f the s tar. A similar result was achieved 
by iGautschvl (|2009l) which is in very good agree ment 
with the spectroscopic study of lLefever et al.l (|2007[ ) who 
placed their sample of BSGs on the instability strip of 
various classes of pulsating stars and verified that these 
super giants are unstable against g-modes (see also I Said 

3.4. Excitation of Nonradial Gravity Modes by the 
e- Mechanism 

From Eq. 22.1 in lUnno etahl (|1989f ) for the hnearized 
equation of energy conservation we have 

\TdynJ \CpJ Lr \ dMrJ 

dlnT/dlnr T ' r dAf/ ^ ' 
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where lo = an T^yn is the dimensionless eigenfrequency, 
dS is the Lagrangian entropy perturbation, is the tan- 
gential eigendisplacement, T' is the Eulerian tempera- 
ture perturbation, and all other symbols have their com- 
mon meaning. In the dense interior of the model, the 
ratio of the thermal time scale Tth oc MrCpT / Lr to the 
dynamical time scale Tdyn ~ {B? /GMY/'^ is large. Even 
for small w, the right hand side of Eq. [7] has to be large 
enough to overcome the high adiabaticity in the interior. 
This expectation is fulfilled for those modes which are 
trapped in the HBSli, for which the perturbation of the 
nuclear energy generation Jcnuc is high enough, and other 
non-adiabatic terms can be overwhelmed. 

To illustrate how the e-mechanism is capable of ex- 
citing g-modes. Figure |4] shows two selected solutions, 
an octupole gov (left) and a dipolc gst (right). The for- 
mer is overstable r] — -1-0.630 while the latter is stable 
77 = —0.208. The radial component of the eingendis- 
placcment £,r{'r) which is shown on top panels (a) is set 
to unity at the surface ^r{R) = 1-0. The remarkable 
feature of gov is its relatively high amplitude in the ra- 
diative zone below the ICZ. There is a phase lag between 
the real and imaginary parts of (r) which accounts for 
the huge peak in the modulus of ^r{r) as seen in Fig- 
uremia. The energy provided by the e-mechanism when 
combined with large amplitude of gov in HBSh drives the 
mode and helps it overcome the radiative damping. In 
panel (b) , the run of energy- related quantities are plotted 
in arbitrary units. The abscissa is logT. The solid line 
is the cumulative work integral as in Eq. [5l the dotted 
line is the work derivative dW/ d In r and the dashed line 
is the square of the Lagrangian entropy perturbation SS 

SS^C,{^^-vJj). (8) 

After a sharp rise in the HBSh, the profile of W flattens 
until it reaches the surface. Above the base of the ICZ, 
the work integral and its derivative are constant, and 
there is no heat exchange in the envelope, so there SS ~ 
0, for this specific solution. Panel (c) shows a zoomed-in 
view of the HBSh capable of exciting g-modes. Clearly, 
the heat exchange and work derivative are in phase. The 
maximum of (SS)'^ and work derivative coincide with the 
location of HBSh where on top of all ^(r) is large, and 
it is the combined constructive effect of all these facts 
which destabilizes gov. The e— folding time Tg for gov and 
other excited modes are less than a year. This means 
that given an initial trigger, pulsational instabilities can 
develop fast (compared to the evolutionary time scales) 
and saturate to reach an observable level on the surface. 

On the right panel, gst has large enough amplitude in 
the HBSh to provide enough energetics for dcstabiliza- 
tion, but the radiative damping above the ICZ dissipates 
the mode energy so that the net cumulative work is neg- 
ative, and this mode is stable 77 = —0.208. Note differ- 
ent values of eigendisplacement between the right and 
left panels. The Fe-bump of opacity embodies negligi- 
ble amount of mass (see Section I2.2p and cannot inject 
enough energy to destabilize gst- 

A noticeable difference between the aforementioned g- 
modes in Figure U] is that the net kinetic energy Ek of 
the overstable mode is twelve orders of magnitude larger 
than that of the stable mode. The simple explanation 
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Fig. 5. — An example of a quadrupole g-mode wrhich is stabilized 
by ft-mechanism. The frequency and growth rate of this mode is 
selected from the broad peak in Figure [3] The Iron-bump has a 
slight stabilizing effect, but the He-bump contributes positively to 
the mode driving; the net effect, however, is mode stability. 

for this is the larger amplitude of gov at HBSh compared 
to gst (see Eq. 2]). That is why those regions in the 
model that are capable of damping gov have no effect on 
its net work W, while the radiative zones have enough 
energetics to dissipate the mechanical energy of gst and 
damp it out. 

3.5. Any Hope on k— Mechanism? 

The broad maxima of stable modes in Figure [3] 
correspond to gravity modes stabilized by the k- 
mechanism arising from the bump of opacity associated 
with the He partial ionization zone at logT w 4.6. This 
zone is situated in the outer dilute envelope of Rigel and 
contains negligible mass (see the long-dashed curve in 
Figure [1]). As a result, the He partial ionization zone is 
incapable of attaining adequate inertia to drive gravity 
modes. 

Figure [S] illustrates an ^ = 2 g-mode selected from the 
broad maxima of stable modes depicted in Figure [3] with 
the frequency of ^0.26 d~^ and 77 = —0.154. The profile 

(,rir) shown in Figure [5] panel (a) has a relatively small 
amplitude in the radiative zone between the Iron and He 
opacity bumps. This is the reason for the negative work 
derivative dW/d In r inward of the He bump (dotted line 
in Figure [5] panel b). However, there is still a positive 
and promising contribution to the work integral arising 
from the He bump as shown by the dotted line in panel 
(b). 

4. CONCLUSION 

The hydrogen burning shell in models of blue super- 
giant stars (like Rigel) is partly situated in the the ra- 
diative layer below the intermediate convective zone. 
This is providing a situation where some modes can 
have a chance to overcome the radiative damping by 
the e-mechanism provided they have large amplitudes 
(or equivalently high kinetic energy Ek) in the radia- 
tive part of the hydrogen burning shell. This can occur 
when the heat exchange (perturbation to entropy change 
(SS)'^) reaches a maximum value at the H burning shell. 
When the location of maximum heat exchange between 
different stellar layers coincides spatially with the loca- 
tion of appropriate fluctuations of the generated nuclear 
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energy, an instability is likely to take place. Our cur- 
rent study which includes effects of rotational mixing 
can explain the long-term (21 ^ P < 127 d) period- 
icities in Rigel. The short-period limit (21 day) where 
e-mechanism destabilizes some g-modes can possibly be 
altered by improvements in the physics of rotating stars 
such as the efficiency of rotational mixing or the extent 
of convective overshooting from convective zones. 

Furthermore, the close similarity between Rigel as a 
late B supergiant and other early A type supergiants 
in photometric and spectroscopic variability suggests ex- 
trapolating the application of the e-mechanism to the 
long-term variation of the a Cyg class of pulsating stars. 

The observed short-term spectroscopic variations of 
Rigel with periods below ^10 days are still theoretically 
challenging (see also iGa utschv 2009.) . Previous works of 
iDziembowski fc Pam iatnvkh (19931): Ibziembowski et all 
(|1993f l and lPamvatnvkh (1999 ) showed that the Iron en- 
hancement and the resulting increase in the height of 
opacity peak can help explain the overstability of main 
sequence pulsating B stars against heat driven p- and g- 
modes. If the observed variations with periods shorter 
than about a week originate from Rigel's pulsations, not 
from other sources (e.g. spots, variable winds, and prop- 
agating shocks), then At-mechanism is the only plausible 
means to induce instabilities. Further local enhancement 
of Iron-group elements at logT w 5.2, in addition to He 
enhancement at logT « 4.6 can possibly lead to exci- 
tation of shorter period modes. However, a more rigor- 
ous theoretical investigation of possible mechanisms to 
locally enhance the desired elements is called for. 

With the advent and success of the current space-based 
missions such as MOST, CoRoT and Kepler, (almost) 
uninterrupted high precision photometry over long time 
baselines is made possible. Supergiants such as Rigel and 
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a Cyg type variable stars with possible long periods will 
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wavelength photometry as plan ned with BRITE nano- 
satellite constellation mission (jKuschnig et al.l [2009). 
The first of these nano-satellites is expected to be 
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pergiants (as progenitors of core collapse supernovae and 
their cores as seeds of the future neutron stars or black 
holes) to come to fruition, and bring a wealth of infor- 
mation about the near-core conditions of the pre collapse 
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